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ABSTRACT 

An acoustic energy reflection experiment was conducted in a labora- 
tory tank containing water and an artificial kaolin sediment. The 
reflection coefficient for the sediment was measured for a limited 
range of angles of incidence at one frequency and one set of measured 
sediment mass physical and viscoelastic properties. A mathematical 
reflection model was derived in terms of the kno^m sediment properties. 
The measured and the theoretically calculated reflection coefficients 
are compared. The measured reflection coefficient and a probable angle 
of intromission showed poor agreement with theoretically predicted 
values. For the sediment parameters of the relatively soft clay used 
here, the reflection model is quite sensitive to values of sediment 
density and sound speed and is only mildly sensitive to values of 
dynamic rigidity and sound absorption coefficient near the angle of 
intromission . 
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I . INTRODUCTION 



A. BACKGROUND 

The development and use of bottom bounce sonar and long range 
passive detection systems, and the use of sonar in shallow water all 
require knowledge of the effect the ocean bottom has on sound propa- 
gation loss. This effect is most often characterized by the quantity 
bottom loss (BL) in decibels, which is related to the reflection co- 
efficient (R) , the ratio of the reflected acoustic intensity to the 
incident acoustic intensity, by the equation BL = -10 log^^ R, The 
large variations in bottom loss throughout the world ocean bottoms, 
coupled with the passive or active sonar figure of merit sensitivity 
to small changes in this parameter, dictate a need for accurate 
measurement of bottom loss characteristics in each area of the oceans 
where bottom sediment parameters are different, 

A direct method for collecting such information would be an 
oceanographic survey ^ of massive proportions, to determine the acoustic 
reflection coefficient at various angles of incidence and various 
frequencies of interest in each distinct bottom region of the world's 
oceans. On the other hand, surveys of other physical characteristics 
of the ocean bottoms have been going on for years, with the result 
that some core samples and raw and analyzed data exist for practically 
every oceanic bottom region. Although much information is available 
in the literature concerning the relation of reflection and transmission 
of acoustic energy to the properties of the reflecting and transmitting 



7 



media, a round table of exploration geophysicists concluded recently 
that their understanding of the "nature of seismic reflection is far 
from complete."* In anti-submarine and submarine warfare, where the 
exact nature of the reflection and transmission of acoustic energy 
at the ocean bottom must be v/ell known, models for calculation of re- 
flection coefficients generally are limited to consideration of a plane, 
homogeneous, compressional wave, transmitted at various angles of in- 
cidence in water and reflected from and refracted in layered sediments 
having various densities and other physical properties, including 
acoustic energy attenuation. 

The energy loss a plane compressional wave undergoes upon reflection 
from the ocean’s bottom is examined in a number of reports, papers and 
textbooks. Three of these references were chosen to describe the nature 
of the problem examined in this paper. White [1] was used as a basic 
text which describes the energy loss as a result of both compressional 
and shear waves being refracted in the sediment. Hamilton [2] uses the 
physical properties of density, compressional and shear wave speeds, 
along with related values of bulk modulus (K) and modulus of dynamic 
rigidity (y) in the sediment to predict values of reflection coefficient 
and bottom loss. Hamilton treats the sediment as both an elastic and 
a viscoelastic medium. Bucker et al [3] compares measured and theoretical- 
ly calculated values of bottom loss for multi-layered viscoelastic 
ocean sediments. 



* 

Anstey, N. A., and Allen, S. J. , "What Direction Should We Set 
for Hydrocarbon Exploration?" Exploration Geophysics , Vol. 21, No. 3, 
p. 412-421, September 1973 
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A number of the acoustic reflection experiments that have been 



conducted at sea, such as those referenced in [2] and [3] , have shown 
varying degrees of agreement between theory and experiment. The dis- 
agreement usually is greatest near the critical angle or angle of in- 
tromission. The theoretical modeling of these experiments has been 
complicated by sediment layering and other non-uniformities and by 
difficulty in obtaining realistic physical properties. A laboratory 
experiment involving a smooth-surfaced homogeneous sediment and an 
undisturbed water layer helps to eliminate some of these complications 
in comparing measured and theoretical reflection coefficients. 

B. DYNAMIC RIGIDITY 

A key element in improved acoustic reflection models has been the 
inclusion of the dynamic rigidity and dissipation properties of the 
sediment. Most researchers, such as Hamilton [ 2 ] and Bucker et al [3] , 
used measured Stoneley wave speeds to calculate the real part of dynamic 
rigidity. Only estimates of the imaginary part were available. Also, 
Stoke *s assumption that in a viscoelastic fluid tlie bulk viscosity is 
equal to zero was employed by these researchers to estimate a value of 
the imaginary part of the first Lame constant (X). 

At the Naval Postgraduate School, the torsional oscillator method 
has been used to measure the real and imaginary parts of the complex 
dynamic rigidity. Cohen [4] described one set of these measurements 
for a saturated, water-kaolin artificial sediment, and Martinek [5] , 
conducting research on the same type of sediment, measured values of 
absorption coefficient. 
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C. REFLECTION COEFFICIENT 

An objective of this research was to examine the usefulness and 
validity of the values of dynamic rigidity and absorption coefficient, 
along with other data on the physical properties of the sediment, in 
calculations with a viscoelastic reflection model, and to compare the 
calculated results with experimental reflection coefficients. The 
viscoelastic reflection model of Bucker et al [3] was modified to avoid 
the use of the assumption that the bulk viscosity is zero for a compres- 
sional wave. Some preliminary experimental measurements were made of 
the reflection coefficient at one frequency and for one sediment condi- 
tion for comparison purposes. 
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II. THEORETICAL REFLECTION COEFFICIENT 



A. DESCRIPTION OF THE MATHEMATICAL MODEL 

Figure 1 illustrates a plane wave reflection model which was 
assumed to be an approximation to the experiment. The symbols and 
notation defined in Figure 1 were drawn from Bucker et al [3] . Refer- 
ring to this figure, (() and m are potential functions defined in the 
following paragraph and x is positive to the right while z is positive 
downward. 

A summary of the development of water-sediment interface equations 
in [3] follows. The wave equation is given by the relationship 
V^(f) + ^ (p = 0, and its assumed potential function solutions are; (1) 

^*2 ” ^*2 ^ 

^*1 ” ^*1 i (kx + a^^kz - tot), (3) 

^”2 ^"2 “ ^2^^ ~ 

y* = B» exp i(kx 4* b kz - tot) (5) 

J- 1 1 / 

= 0 , B"^ = 0 . 

In water only the real values of k are considered for 

(0 sin 0. 

k = i . (6) 

c 

w 

The term c^ is the measured speed of a compressional wave in water, 

A * 2 A ”2 are the measured amplitudes of the incident and reflected 

waves, respectively, A*^ and B*^ are unknown quantities, 0^ is the angle 
of incidence, and c[)"^ and are both equal zero. 
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Figure 1. Theoretical Reflection Model. 



For the v;ater-sediment interface two boundary conditions may be 
applied to the potential function: 

1. the vertical particle velocities (w) are continuous across the 
interface; 
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(7) 



2. the vertical stress- tensor elements are continuous across the 
interface; 



i / 2i 

= 77(-PW - 


a^(}) , 

2y — - + 2]i ) 


(8) 


ZZ 0) 


9x^ 9x9^ 




p = 0 . 






zx 






Following Bucker [3] , three interface equations are derived as 


follows 


from the boundary conditions and 


treatment of the sediment as 


a visco- 


elastic half-space; 






' a^(A’ - A" ) 


= a^ A* + B’ 


(9) 


2 2 2 


11 1 


P2 <^'2 " ^" 2 ’ 


= y^[(b^^-l) (A'^) +2b^(B'^)] 


(10) 


2a^ (A'^) 


= (b^^-1) B'^ - 


(11) 


where: 


2 




2 


Pi c 


(12) 


"1 


2y^ " ^ 




2 




b/ 


- -1 


(13) 


1 





13 



c 

w 

c = — : = horizontal phase speed in water 

sin 0. 

1 

= cot 0. 

2 1 

= density of sediment 
= density of water 



The Lame constants 


A and p are comi^lex, 




II 

\ 


iA" 


(14) 


1 1 


1 




b = - 


«— 1 
•H 


(15) 



B. SEDIMENT PARAMETERS 

The above model permits a solution to the wave equation in terms of 
sediment parameters A, y, p and c, of v/hich A is the only term which 
does not have independently determined values. Cohen [4] measured 
values of dynamic rigidity, y*^ and y”^ for a range of sediment densi- 
ties. Martinek [5] measured the absorption coefficient, a^, and the 
compressional wave speed in the sediment, c^^, at a sediment density 
comparable to Cohen’s and at frequencies in the range of experimental 
interest. These values are summarized in Table I of Section IV. Addi- 
tional sound speed and density values were measured as part of this 
experiment and recorded in Table II of Section IV. 

In the absence of scattering from the water-sediment interface, a^ 
in theory is related to A"^ and y*'^ and c^ is related to A*^ and 
With these known quantities, the next series of steps transposes the 
equations to eliminate A * and A"^ as independent variables and solves 
for a reflection coefficient in terms of sediment parameters y*j^/ ^”l' 
Pf/ along with c^, p^ and 0^ in the water layer. 



14 



/ 







C. CALCULATION OF REFLECTION COEFFICIENT 

In the following development, the symbols are similar to those used 
by Borcherdt [6] but with a reversed sign convention for the complex 
exponential. The sign convention used here is the same used by Bucker 
et al [3]. Definition of terms is as follows; 



k = k* - ik" ^ = complex propagation constant 
a = a* - ia" = complex wave speed and 



00 = angular frequency 



a = p.v 




where is the bulk modulus and p.v. indicates principal value 
From Hamilton [2] , ~ that 




1/2 



(16) 



The compressional v;ave speed in the sediment, c^, is given by 



c 



1 



1 



= Re [i] 



(17) 



The compressional wave absorption coefficient in the sediment 



IS 



a = oolm [ — ] 
p a 



(18) 



The complex compressional v/ave propagation constant is 




It 
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For physical reasons, > 0 to ensure decreasing values of the 
potential function with positive increasing values of distance along 
the propagation path. By the same physical reasoning, > 0. 

Algebraic steps to solve for reflection coefficient, R, in terms of 
y*l, frequency, f, and independent of and are 

shown in detail in Appendix A. These steps yield the following set 
of equations; 



R E 





A-2 


2 


Z-1 








Z+1 



Z = 



^'2 ^"2 

A*2 - A-2 



2a^y^k 



p2(bi^+l)W^ 






(b ^-1)^ 



/ \ 



a^ = H + iJ 



= F + iG 



F = 



G = 



2 , 

P c u 

Pf M 1 






- 1 



2 „ 



(y'^)^+(y"^)^ 



(A-10) (A-ll) 



(A-9) 

(A-12) 

(A-15) 

(A-16) 

(A-17) 



4'.f 



- I -1 



J = 



TT.f 



c a 2 
1_ (—-£-) 
^2TTf ’ 



H + 1 



(A-18) 



(A-19) 
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Equation (A-11) for R then was solved using standard WATFORG on 
the IBM 360 computer at the Naval Postgraduate School. Values of 0^ 
were chosen from zero to 7T/2 and the physical properties of the 
sediment were taken from Cohen [4] and Martinek [5] . Values for BL 
and R were computed over the range of experimental error measured in 
the physical properties cited in these two references, in order to 
examine them for sensitivity to changes in these values. 
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III. MEASUREI'OENT OF REFLECTION COEFFICIENT 



A. SUr^RY OF THE EXPERIMENT 

The experiment was conducted by reflecting a pulsed beam of sound 
from a water-sediment interface at various angles of incidence. The 
incident and reflected wave amplitudes were detected with a directional 
receiver of design identical to that of the transmitter. A sound fre- 
quency of 120 kHz was used. The received signals were cimplified, fil- 
tered and signal processed to improve signal-to~noise ratio. The 
reflection coefficient v;as computed from the ratio of the squares of 
the voltage amplitude of the direct and sediment reflected pulses, 
using appropriate corrections for spherical divergence. 

B. THE APPARATUS 

The artificial sediment, a commercially available pure kaolinite, 
was contained in a large, rectangular tank, 2.42 m long by 1.20 m wide 
by 1.00 m deep. The leveled, saturated sediment thickness was 60 cm. 
The water depth was 31.5 cm. 

Identical transmitting and receiving transducers were set up in 
the water and mounted on a horizontal rail above the tank. Each trans- 
ducer could be rotated through 360° independently in both the vertical 
and horizontal axes and moved horizontally and vertically in the tank. 
The two electrostatic type transducers, built by research physicist 
Donald F. Spiel of the Naval Postgraduate School, have an active circu- 
lar face of 6.0 cm diameter, consisting of an aluminum coated, one mil 
thick plastic film mounted over an aluminum electrode which is set in 
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a larger plastic disc. In order to keep the air gap between the plastic 
film and the electrode uniform, equalization of pressure betv;een this gap 
and an air filled space behind the electrode was achieved by small holes 
drilled in the aluminum electrode. The inside of the transducer was 
partially evacuated. At 120 kHz the 3 dB beam width is about 20®. 

The basic geometry of the experiment is shown in Figure 2. Referring 
to this figure, both the transmitting and receiving transducers were 
rotated to a selected angle of incidence, 0^. 

Based on the choice of centerline height of each transducer being 
23.5 cm above the sediment layer, the idealized ray path y^ was computed 
for each 0^^ and the transducers were positioned along the horizontal 
rail so that y^ = y^ and “ ^ 2 * horizontal x^lacement also was 

chosen to ensure that the same region of sediment v/as ensonified for each 
case. 

The range of angles accessible for this exiDeriment is limited by 
the dimensions of the apx^aratus and the necessity to keep the sediment 
reflecting surface in the far field region of the directional trans- 
ducers. The extent of the far field region was determined experi- 
mentally by plotting the directly propagated sound field pressure as 
a function of range. 

The directivity of both transducers at 120 kHz provided a suffi- 
cient reduction in off-axis signals to eliminate interference from 
side wall and water-air surface reflections. This was confirmed experi- 
mentally. At angles above 70® a small amount of interference between 
direct and sediment reflected pulses could be observed, the effects of 
which are discussed later. 



19 




Figure 2. Experiment Geometry. 
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C • INSTRUMENTATION 



Figure 3 illustrates the transmitting and receiving equipment and 
circuitry in block diagram form. This configuration is similar to 
that used by Martinek [5] . The Hewlett Packard 204C oscillator fre- 
quency was set at 120 kHz within + 3 Hz and was monitored with a fre- 
quency counter. The General Radio 1396-A tone burst generator was 
set to produce two-cycle tone burst pulses with an approximate 10 msec 
delay between pulses. This delay v;as sufficient to prevent undesired 
reflections from disturbing the measurements. The Hewlett Packard 
467 a amplifier was used to boost the voltage by a factor of ten before 
the signal went to a 1 :3 transformer for further amplification. The 
polarizing circuits for both the electrostatic transducers are shown 
in Figure 4. 

The received signal was amplified either 20 or 40 dB by a Hewlett 
Packard 465A amplifier and then passed through Krohn-Hite high and low 
pass filters set to a 20 kHz pass band centered at 120 kHz. A Prince- 
ton Applied Research (PAR) model 160 boxcar integrator was used for 
signal processing. The time delay feature of the boxcar integrator 
was used to determine the times of pulse arrival. The enhanced signal 
from the arriving pulse was then fed to a Moseley X-Y Recorder. Fig- 
ure 5 which was traced from original records represents typical received 
wave forms after integration and recording o Wave form points selected 
to determine zero to peak voltages are noted on this figure* 

D. DATA COLLECTION PROCESS AND ERROR ANALYSIS 

A number of angles of incidence were selected and appropriate hori- 
zontal separations were pre-calculated. At each setting of horizontal 
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Figure 3« Instrumentation. 
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Figure 4« Polarizing Circuits. 
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Figure 5. Comparison of Direct Path and Reflected Wave Forms. 



I 

I 

spacing, measurements were recorded of the received pulses for both a 

direct path between source and receiver (6. =0) and for the sediment 

1 

reflected path. Additional gain was required for the latter measure- 
ment. The receiver then was moved its full horizontal travel to 
determine that no other pulses were being reflected from belov; the 
water-sediment interface . 

Referring to Figure 5, zero to peak voltage was determined by 
averaging the maxiraiim positive and negative peaks of each recorded 
pulse. In all cases these peaks occurred in the third cycle of the 
pulse. Voltage output anomalies in the tone burst generator, trans- 
formers, filters, amplifiers, boxcar integrator and the X-Y recorder 
were assumed to be the same for both direct path and reflected path 
voltages and therefore were cancelled by division in the reflection 
coefficient ratio. Also, because of this ratio computation, trans- 
ducer sensitivities, sound pressure level calculations, etc., were 
not required. 

The log-log plot of direct path average zero to peak amplitude ver- 
sus horizontal separation v;as constructed to determine the minimum 
horizontal separation at which the pressure varies linearly with the 
inverse range (1/r) . This established the inner limits of horizontal 
separation which in turn set the minimum angle of incidence at 
^i ~ • Since the reflected path was longer than the direct path, 

this separation insured that the ensonified patch of sediment was also 
in the far field. The dimensions of the tank limited the maximum 
angle of incidence to 0^ = 76°. 

The arrival time of the sediment reflected pulse occurred later 
than the arrival time of the direct pulse. The measured arrival time 
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differences between the two pulses were compared with time differences 
calculated using the idealized geometry ray path difference and the 
known speed of sound. Compatibility between observed and calculated 
delays confirmed detection of the sediment reflected signal. 

The uncertainty in the angle of incidence (0^^) on each transducer 
is about ^ 1.0°. Part of this is due to the precision with which the 
scale could be read and part arises in the assembly process of aligning 
the transducer face to the zero mark on the angle indicator dial. It 
is estimated that horizontal separations were accurate to + 0.2 cm. 

The height of the transducer above the ensonified sediment patch was 
accurate to +_ 0.2 cm and the sediment surface within the ensonified 
patch v/as level and smooth within ^0.1 cm. With a sound wavelength of 
about 1.2 cm, this scale of roughness in bottom level is considered not 
to have a significant effect on observed signal amplitudes. The un- 
certainty in measuring voltage amplitude is estimated to be 2% of the 
millivolt scale in Figure 5. 

At angles of incidence greater than 70°, small but perceptible 
off-axis energy from the direct path was detected. Since the 
difference between reflected and direct path lengths decreases with 
greater angle of incidence, the arrival time difference between the two 
pulses decreased so that the leading edge of the reflected pulse over- 
lapped the tail of the direct path off-axis pulse. This resulted in 
an estimated + 5% error in reflected voltage £implitudes above 70°. 

E. THEORY OF MEASUREMENT 

Referring to Figure 2, it was assumed that the sound pressure ampli- 
tude varied inversely as the square of the range from the source and 
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that the absorption in the water was negligible. The assumption was 
also made that the sediment surface was sufficiently smooth and homo- 
geneous so that the reflected rays behaved as if they had originated 
at the mirror image position of the source. 

The received intensity I, at any range x along the direct path, was 
determined by orienting the beam axes of the transducers toward each 

other and measuring signal amplitude A*. Intensity I is proportional to 
2 

(A*) . The reflected wave intensity which is proportional to the 

square of the measured reflected wave amp] itude A*, is a function of 

the reflection coefficient R, and the spherical spreading along 

the total reflected path length y = y^ + y2 = 2y^, assuming that the 

reflection is specular. is the axial intensity at unit distance 

from the sound projectorc Stated in equation form, this relation is 
2 2 

= (A") = Rl^/y . Solving for R and substituting for determined 

from the direct path measurements gives the following: 




Bottom loss is calculated by 



BL = -10 log^^R (dB) 
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IV. PRESENTATION OF DATA 



Tables I and II, which list some of the mass physical prox:>erties 

of the sediment, are self-explanatory. Table II records the speed of 

sound and density in the sediment for three different dei^ths near the 

surface. Sound speed ratio o^/o was measured v/ith a sound velocimeter 

1 w 

and densities were measured using the sampling device described by 
Martinek [5] . The data show both density and velocity gradients in 
the upper part of the sediment. 

Table III lists a series of eight combinations of the five measured 
sediment parameters which are used in the tlieoretical model. These 
combinations are designed to determine the model's sensitivity to para- 
meter changes with each combination given a sequence number for refer- 
ence purposes. Sequence (1) is based on Table I dynamic rigidities at 
a density of 1.35 gm/cm^ along with sound speed and absorption coeffi- 
cient at a density of 1.38 gm/cm^. The ma>:imum positive experimental 
error from Table I dynamic rigidities is the cliange inserted in Sequence 

(2) . The effect of only a density change from 1.38 gm/cm^ to 1.28 gm/ 
cm^ is determined from using Sequence (3). Together, Sequence (2) and 

(3) test the model for sensitivity to density changes. Sequence (4) 

3 

shifts to the Table I values of dynamic rigidity for 1.28 gm/cm density 

and in Sequence (5) the measured sound speed for a density of 1.28 gm/ 

3 

cm in Table III is inserted. Together, Sequence (4) and (5) determine 
the effect of only a sediment sound speed change. Also, Sequence (5) 
is based on the measured data which is assumed to be the closest approxi- 
mation to the sediment conditions during the reflection experiment. 



28 



I 



Sequence (6) and (7) are designed to test the model for large arbi- 
trary changes in dynamic rigidity. The density of 1.24 gm/cm^ is an 
assumed value of sediment density at a point closer to the surface 
than that at which a suitable sample could be collected. Sequence (8) 
tests the model for a large arbitrary change in absorption coefficient 
relative to that in Sequence (5) . 

Tables IVa and IVb tabulate the values of R and BL for the eight 
sequences computed from Equation (A-11) for a range of angles of inci- 
dence from 14® to 83®. From these two tables sensitivity of the 
theoretical model is determined. The values from Sequences (1) and (5) 
are plotted in Figure 6 for comparison. 

Table V tabulates the data from the reflection experiment. The 
values of BL in this table are plotted in Figure 6 along with calcu- 
lated values from Sequences (1) and (5) . 



29 



I 

i 



TABLE I 



Physical Properties of Saturated Water-Kaolin Sediment 



Cohen [4] Sediment Properties for 20®C 



Density 


Dynamic Rigidity 


Pi 


^'i 


l."l 




^-4 


“4 




X 10 


X 10 


(gm/cm^) 


2 

(dynes/cm ) 


2 

(dynes/cm ) 


1.28 


26.6 + 3.4 


2.85 + 1.20 


1.31 


54.8 + 7.4 


9.68 + 4.36 


1.35 


67.9 +14.5 


16.6 6.90 









Martinek [5] Sediment Properties for 19,9®C 



Frequency 


Density 


Sound Speed 
** 


TVbsorption 

Coefficient 


f 


p. 


c. 


a 




■^1 


1 


P 


(kHz) 


3 

(gm/cm ) 


{ cm/sec) 


(nepers/cm) 






2 


-3 


100 


1.38 


1422 X 10 


3.79 X 10 






2 


“3 


125 


' 1.38 


1422 X 10 


4.83 X 10 






2 


~3 


* 120 


1.38 


1422 X 10 


4.62 X 10 



Linear interpolation from 100 and 12b kHz data 



c^/c = 0.94 
1 w 
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TABLE II 



Measured Sound Speed and Densities 



Average 

Sediment 

Depth 


Average 

Sediment 

Density 


— at Depth 
c 

w 


** 


(cm) 


(gm/cm^) 


(cm) 


(cm/sec) 


1.55 

(Surface sample) 


1.28 + .05 


0.974 at 1.37 


1442 X 10^ 


3.85 


1.38 + .05 


0.965 at 3.39 


1428 X 10^ 


7.85 


1.42 + .05 


Not measured 


- 



* 

Measured at 18.8°C 



* * 2 

Corrected to 20°C, based on c = 1480 x 10 cm/sec 



TABLE III 



Combinations of Sediment Parameters 
Used in Theoretical Computations 



Sequence 

No. 


y, 1 

(dynes/ 


cm ) 


'’l 3 

(gm/cm ) 


(cm/sec) 


a 

P 

(nepers/cm) 


1. 


67.9x10^ 


16. 6x10"^ 


1.38 


Io422xl0^ 


0.462x10“'^ 


2. 


82.4xlo'^ 


23.5x10"^ 


1.38 


1.422x10^ 


0.462x10“'^ 


3. 


82.4x10^ 


23.5x10^^ 


1.28 


1.422x10^ 


0.462x10“"^ 


4. 


26.6x10“^ 


2.85x10"^ 


1.28 


1.422x10^ 


0.462x10“^^ 


5. 


26.6x10^^ 


2.85xl0'^ 


1.28 


1.442x10^ 


0.462x10“'^ 


6. 


266 


28.5 


1.24 


1.450x10^ 


0.462x10"'^ 


7. 


26.6x10^ 


2.85x10^ 


1.24 


1.450x10^ 


0.462x10""^ 


8. 


26.6x10^ 


2.85xl0'^ 


1.28 


1.442x10^ 


0.462xl0"^ 



31 



n 






I 



Calculation of Reflection Coefficient and Bottom Loss 

Sequences 1-4 
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Sequences 5-8 
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Measured Reflection Coefficient and Bottom Loss 



I 



Bottom 

Loss 

BL 

(dB) 


41.5 

42.0 

42.9 

43.0 

42.0 
42.3 

40.9 
38.2 


Reflection 

Coefficient 

4 

R X 10 


0.715 
0.639 
0.505 
0.500 
0.628 
0.578 
. 0.822 
1.506 


Reflected 
Path 
0-Peak 
Avg. volts 

(mv) 


0.625 

0.550 

0.440 

0.375 

0.335 

0.300 

0.305 

0.375 


Direct 
Path 
0-Peak 
Avg. volts 

(mv) 


96.5 

84.0 

71.5 

58.5 

45.0 

41.5 

35.0 

31.5 


Path Length 
Difference 

(cm) 


17.10 
■ 14.82 

12.60 
10.82 
8.28 
7.44 
6.62 
5.78 


Reflected 
Path Length 

Yl+Yj 

(cm) 


73.12 

81.94 

94.00 

111.62 

137.42 

152.10 

170.52 

194.28 


Horizontal 

Separation 

X 

(cm) 


56.02 

67.12 

81.40 

100.80 

129.14 

144.16 

163.90 

188.50 


Angle of 
Incidence 

e. 

1 

(degrees) 


50 

55 

60 

65 

70 

72 

74 

76 
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Bottom Loss (dB) 
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V. COI^ ARISON OF DATA 



A comparison of the bottom loss calculations from the reflection 
model with the measured results from the experiment reveals somewhat 
poor agreement between the two. The experimental limits on the angles 
of incidence prevented measuring reflection* coefficients in the 
vicinity of the theoretically predicted values of 40° to 43°, Figure 6 
shows a possible measured angle of intromission in the vicinity of 65°, 
but the lack of experimental results below 50° leaves open the possi- 
bility of even higher measured bottom losses at the lower angles. An 
additional comparison was made with a liquid-liquid model described in 
Kinsler and Frey [7, p. 145], using only the densities and sound speeds 
of Sequences (1) and (5) , which predicted an angle of intromission in 
the range of 26° to 27°. The theoretical magnitudes of bottom loss 
near the angle of intromission are 3 0 to 20 dB higher than the values 
measured at about 65°. 

The eight combinations of parameter changes reveal a number of 
facts regarding the theoretical model and its application to the experi- 
ment performed. A summary of comparisons among sequences is contained 
in Table VI on the following page. The high sensitivity to density and 
sound speed changes shown in this table agrees with what is expected 
from the basic liquid-liquid model [7] . The Table VI comparisons also 
reveal that dynamic rigidity and absorption coefficient affect bottom 
loss magnitude only near the angle of intromission. The range of para- 
meter changes selected for the sensitivity tests had only a small effect 
on the range of theoretical angles of intromission, as shown in Table VI. 
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TABLE VI 



Comparison of Sequences 



Sequence 


Parameter 


Sensitivity of 


Numbers 


Change 


Equation (A-ll) 


1 and 2 


Small change in 
dynamic rigidity 


Sensitive only at 
the angle of intro- 
mission 


2 and 3 


Density 


Very sensitive 


4 and 5 


Sound Speed 


Very sensitive 


6 and 7 


Large change in 
dynamic rigidity 


Sensitive only near 
the angle of intro- 
mission 


5 and 8 


Absorption 

Coefficient 


Low sensitivity only 
at angle of intromis 
sion 



37 



VI. CONCLUSIONS AND RECOMMENDATIONS 



The data collected for one sediment condition at one frequency for 
a limited range of angles of incidence must be considered as preliminary 
in nature. Three extensions of data collection using the developed 
experimental method are needed to complete •the experiments. A variation 
in sediment conditions is needed to measure the effect on reflection 
coefficient from known changes in density ^ sound speed and dynamic 
rigidity. Apparatus now exists to measure all of these parameters. A 
variation in frequency also is needed to measure the effect of changes 
in absorption coefficient, which has been assumed to be the only fre- 
quency dependent sediment property. Measurements over a wider range of 
angles of incidence are needed in order to determine experimentally the 
existence of an angle of intromission. In order to increase the range 
of angles of incidence changes in the appai'atus will be required. As 
an example, a deeper and larger tank will permit measurements at smaller 
angles of incidence. 

The low sensitivity of R to values of dynamic rigidity could be due 
to the range of dynamic rigidities used for the material not being repre- 
sentative. Because of the fact that the reflection model is quite 
sensitive to sound speed and density, there exists the possibility that 
a part of the discrepancy between calculated and observed reflection 
coefficients is due to inadequate accuracy for measured sediment physi- 
cal properties. 

The sediment densities and velocities measured during this experi- 
ment and earlier by Martinek [5] confirm that gradients in properties 
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do exist in the upper portion of the sediment. AJ. though the gradients 
may have a significant influence on the observed reflection, this was 
not accounted for in the theoretical model used here and time limita- 
tions did not permit the adaptation of another model. One such model 
which does include sound velocity gradients in the reflecting medium 
has been developed by Morris [8] based on the theory of layered media 
reflections described by Brekhovskikh [9] . .The gradient models must be 
included for experiments in which the sound reflection process involves 
greater x^enetration of sound into the sediment than was probably the 
case in the exx^eriment rex^orted here. 
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APPENDIX A 



Algebraic steps to solve for reflection coefficient in Section II 
are set forth in tliis appendix. 



,^,2 _ . ,2 ^ 2,1/2 -iO -ie 

(a) = x-iy = (x + y ) e = re 



(A-1) 




Note : 0 not related to 6 . 

1 



X = X' ^ + 2 ji'^ 



y E X"^ + 211 "^ 



- , 2 21/2 
r = (x H- y ) 



(A-2) 

(A-3) 

(A-i) 



, ,1/2 -i 0/2 , ,1/2, 

a = (r) e = (r) (cos 0/2 - i sin 0/2) 



The positive roots of the following trigonometric identities were 
substituted: 



„ ,1 + cos 0/2, 1/2 . _ ,1 - cos 0/2, 1/2 

cos 0/2 = + ( — '—) , sin 0/2 = + ( ^ —) 



a = {[(i) (1+ - i[(i) (1- 

2 r 2 r 



1 r , ,1/2 . , , 1/2, 

[(r+x) - i(r-x) ] 



l/^ 



Rationalizing the inverse of the above to obtain ^ yields the 



following : 



. , X 1/2 . , , 1/2 

1 (r+x) + i(r-x) 



a 



\/2 
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Substituting the real and imaginary parts into Equations (17) 



and (18) yields 



the follov;ing: 



Re[-] = 

a 



1/2 






a 



1/2 



l/2 



C 



1 



r /2 
y/r + X 



a 

P 



y/2 7Tf y/r-X 
r 



(A- 5) 



(A-6) 



Algebraic manipulation of and a then was used to eliminate 

1 P 

r and state all relations in terms of x and y. 



2r 

r+x 



a = 
P 



2(TTf) (r-x) 



2 2 

r + x - ( — -) r = 0 
^1 



r-x - ^(ri ) r = 0 

2 TIT 



2. . - 



r] = 0 



X 

2 

r 



_ i(^) 

4 



(A-7) 
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Manipulating further as follows: 



2 2 
(c,)^(a 
1 p 



(27rf) (r-x) 

r+x 



( J ^)2 

^27Tf ^ 



r-x 

r+x 



X 

r 



c a 

r ^ 



1 + 

2TTf 



From 



r = 



2 2 
X + y , 



- 1 - A' 

r r 






= 1 - 





c, a 


c a 




1-2 


( ^ + 


^TTf ‘ 


4 

1 




c, a 


c a 




1 + 2 


, ] P>2 
27Tf 


( " P) 

^2TTf 


4 

1 







(^) 

r 



= + 



1 + ( 



c a 
1 P 

^f 

c.a 

1 P v2 

27Tf ^ 



and since a^ and are greater than zero, the positive root is chosen. 
Combining x/r and y/r to obtain y/x yields the following: 



y 

X 



c, a 

J- p 

TTf 



1 - 



(A-8) 



Substituting the definitions of x, y and r in Equations (A-2) , 

(A-3) and (A-4) shows that, related only to 

the groupings of known physical properties of Equations (A-7) and (A-8) . 
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Solving the interface e^l^^tions by adding Equation (9) and (1]) 



and substituting the result into Equation (10) gives: 






, (b^ - 1) (b^ -1)B' 

^2 ^" 2 ^ "" ^ll 2a, 



7 



(b^ -1)^ 



Substituting 



B’ = 



Z = 



A* + A”, 



A* — A** 

^2 2 



^"2 <^'2 - ^”2^ 

(b^^ + 1) 



2a2bik 



2 2 

p^W (b ^ + 1) 



2 2 
(b^^ - 1)^ 



2 a _ 



+ 2b. 



(A-9) 



Letting Z equal the right hand side of Equation (A~9) and defining 
intensity reflection coefficient 

A'* 



R = 



AV 



(A-10) 



such that 



A’2 A»^ 1 f A’yA’^ 1 

^'2 ■ ^"2 ^ “ ^"2^^' 2 1 - 



= Z 



then 



R = 



Z-1 

Z+1 



(A-11) 
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I 

I 

I 



I 



h 



All quantities on the right hand side of Equation (A-11) are x^hysi- 
cal properties of the sediment or water and Z is a complex number. 
Formulating a^ and as complex numbers gives the following: 



^ , w , 2 

^sin 0.^ 

1 



1 iA'‘^) + iy"^) 



= H + ij 



(A-12) 



H and J are defined as follov/s: 



H 



f>i 

1 

(A'^ + 2m'^)^+(A"^ + 2\x"^)^ 



- 1 



(A-13) 



J = 



Pi ( 

'^l sin 0. 1 

1 



(A’^ + 2 m'^) +(A"^ 



+ 2m"^) 



(A-14) 



, w ,2 

Pi (ILT0-) 

1 



= F + iG 



(A-15) 



F and G are defined as follows: 



. , w , 2 

pi^i 

F = 2 H 1 

{\i\) + (V\) 



(A-16) 



G = 



fF"i 



A4 



1 



are contained only 



At this point the two Lame constants X' and X" 

1 X 

in the terms H and J. Equations (A-13) and (A-14) are then restated in 
terms of x and y as defined earlier in Equations (A-2) and (A-3) . 



, 2 



P^x {— ^ / 

1 sin 0, p^x c 

H - ■ ^ 1 
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Substituting from Equations (A-7) and (A-8) then yields the final 
expressions for H and J, 
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